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THE  HYDRODYNAMIC  WAKE  OF  A  SURFACE  SHIP: 
THEORETICAL  FOUNDATIONS 


INTRODUCTION 

A  model  for  describing  the  hydrodynamic  wake  of  a  surface  ship  is  developed.  The  model 
accounts,  in  some  manner,  for  all  significant  wake  flow  processes.  Simultaneously,  the  model  is  formu¬ 
lated  to  be  computationally  realizable.  This  is  necessary  because  small  scale  flow  behavior  on  a  large 
spatial  domain  is  characteristic  of  the  wake. 

During  the  course  of  the  derivation,  we  observe  that  several  approximations  must  be  made  to 
obtain  the  model.  The  approximations  result  from  a  lack  of  firm  theoretical  foundations  in  certain 
areas.  Besides  turbulent  flow  dynamics,  these  areas  include  descriptions  of  multiphase  flows, 
turbulence/ free-surface  interactions,  and  ambient  background  characteristics. 

Other  assumptions  and  constraints  are  required  to  achieve  computational  realizability.  These 
include  the  decomposition  of  the  flow  field  into  potential  and  viscous  dominated  parts,  limits  on  allow¬ 
able  ship  motions,  and  certain  linearizations  and  simplifications  of  boundary  conditions. 

The  model  does  not  include  several  flow  phenomena  considered  of  secondary  or  intermittent 
importance.  These  phenomena  are  associated  with  ambient  stratification  and  ambient  surface  films. 

EQUATIONS  OF  MOTION 

Visually,  the  wake  of  a  surface  ship  contains,  among  other  things,  a  foamy  mass  of  white  water 
caused  by  sprays,  turbulent  entrainment  of  air,  and  bubbles.  To  the  modeler,  the  wake  is  a  multiphase 
flow  field  consisting  of  a  gas-water  mixture.  This  leads  to  immediate  modeling  difficulties  because  a 
firm  theoretical  foundation  for  multiphase  flow  dynamics  does  not  exist  [1,2].  As  noted  in  Ref.  1,  this 
results  partially  from  the  particle/continuum  dichotomy  of  multiphase  flows  and  partially  from  the  in¬ 
ability  to  ascertain  fluid  constants  in  multiphase  flow  systems. 

Therefore,  to  make  any  progress  in  describing  the  wake  of  a  surface  ship,  we  must  develop  an 
approximate  treatment  of  the  multiphase  flow  field.  To  this  end,  we  adopt  an  "inertialess  particle" 
hypothesis  [3].  For  our  problem,  the  validity  of  this  hypothesis  requires,  first,  that  the  gas  particles 
(bubbles)  entrained  in  the  water  column  be  dispersed  enough  to  be  essentially  noninteracting,  and 
second,  that  the  inertia  of  the  gas  particles  be  much  less  than  the  inertia  of  the  entraining  water.  In 
ship  wakes,  the  first  requirement  is  observed  empirically  to  be  met  away  from  boundaries.  The  satisfac¬ 
tion  of  the  second  requirement  then  follows  from  the  satisfaction  of  the  first  because  of  the  relative 
density  ratio  between  gas  and  water. 

In  addition  to  the  above  requirements,  the  inertialess  particle  hypothesis  presumes,  from  the  con¬ 
tinuum  viewpoint,  that  any  infinitesimal  volume  of  water  is  still  large  enough  to  allow  a  meaningful 
specification  of  gas  particle  distribution  in  that  volume.  We  denote  this  specification  by  the  gas  particle 
(or  bubble)  volume  distribution  function  f0yB  which  gives  the  number  of  bubbles  per  unit  of  water 
volume,  per  unit  of  bubble  volume  increment  dVg.  Here,  Vg  represents  the  volume  of  a  gas  particle. 
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From  the  definition  of  f0Y  ,  certain  useful  quantities  related  to  the  bubble  population  can  be 

B 

obtained  directly.  In  particular,  the  number  of  bubbles  per  unit  of  water  volume  having  a  gas  volume 
between  VB\  and  VB2  is 


Q2 

VBl  foVBdVB’ 


(1) 


and  the  total  gas  volume  per  unit  of  water  volume,  or  void  fraction,  is 

To b  =  J0  fovB  VBdVB. 


(2) 


Within  the  inertialess  particle  hypothesis,  the  velocity  uqBi  of  the  gas  particle  is  assumed  to  be  the 
velocity  u0i  of  the  entraining  water,  plus  the  rise  velocity  uR3  of  the  gas  particle  in  otherwise  undis¬ 
turbed  water.  Here,  we  have  adopted  indicial  notation  and  defined  the  "3"  direction  as  upwards.  We 
find 

u0  Bi  =  u0i  +  uR$i3  (3a) 

where  8y  is  the  Kronecker  delta.  The  rise  velocity  itself  is  calculated  independently  of  the  water  ve¬ 
locity  by  methods  such  as  those  discussed  in  Ref.  4.  The  position  vector  R0Bi  of  the  gas  particle  is 
obtained  from  the  integral  of 

RoBi  =  uOBi ■  (3  b) 

The  equations  of  motion  for  the  water  in  the  inertialess  particle  limit  are  developed  straightfor¬ 
wardly  by  noting  that  inertialess  implies  that  no  mass  or  momentum  fluxes  are  carried  by  the  gas  parti¬ 
cles.  Hence,  the  primary  influence  of  the  gas  particles  is  to  change  the  density  of  the  water  from  its 
true  value  p  to  an  apparent  value  p(l  —  yoa)-  Substituting  this  apparent  value  for  the  density  in  the 
flux  (or  conservation)  forms  of  the  continuity  and  Navier-Stokes  equations  [5],  we  obtain  the  equations 
of  motion  for  the  water  as 

[(1  ~  yos)p]  4-  —  [(1  —  yoB^Puok^  =  0  (4a) 

at  axoic 


and 


B_ 

dt 


1(1  —  yofiWo/]  + 

axok 


(d  -  yo*)p«o,«oJ 


an0 

3*0/ 


-  (1  -  yoB'tpg&.i  + 


B 

dx0k 


(1  -  y0fl)  P- 


Buy  Buok 
Bx0k  Bx o, 


Bx0j 


(4b) 


In  Eq.  (4b),  n0  denotes  the  total  pressure,  g  the  acceleration  due  to  gravity,  and  p.  the  viscosity  of  the 
water.  The  appearance  of  the  factor  (1  —  y0B)  modifying  p  follows  from  the  inertialess  particle 
hypothesis  [1,  Eq.  (6.16)]. 


BOUNDARY  CONDITIONS 


Consider  the  two  coordinate  systems  shown  in  Fig.  1.  The  00xoiXq2Xo3  system  is  fixed  in  the  fluid 
and  is  the  coordinate  system  in  which  Eqs.  (1)  to  (4)  have  been  derived.  The  undisturbed  free  surface 
coincides  with  the  x03  =  0  plane.  The  Osxslx52%3  system  is  fixed  in  the  ship.  With  respect  to  this 
latter  system,  RSHi  denotes  a  vector  from  Os  to  a  point  on  the  surface  of  the  ship.  For  those  points 
that  have  motions  in  addition  to  the  rigid  body  ship  motions  (for  example,  propellors  and  control  sur¬ 
faces),  Rshi  is  time  dependent. 
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Fig.  1  —  Coordinate  systems  and  notation 


We  also  define  the  vector  from  O0  to  Os  by  R0Si.  Then,  with  respect  to  the  00*01*02*03  coordi¬ 
nate  system,  the  vector  R0Hi  from  O0  to  the  same  point  on  the  surface  of  the  ship  is  obtained  [6], 

Roh,  =  Rost  +  TjjRSHj.  (5) 

The  transformation  matrix  T  is  determined  from  the  noncommutative  set  of  rotations— yaw  0S3  about 
pitch  0S2  about  xS2,  and  roll  0S1  about  x51— with  positive  senses,  as  indicated  in  Fig.  1.  From  Ref. 
6,  we  find 


T  = 


sin  6 si  cos  0S2  sin  053  sin  6S2  sin  0Si  sin  0S3  sin  6S2  cos  i 

+  cos  6 si  cos  —  cos  6 si  sin  6s\ 

—  sin  0S2  COS  052  sin  05!  COS  052  COS  051 


(6) 


The  velocity  Ush,  of  this  point  on  the  surface  of  the  ship  is  given,  with  respect  to  the  ship  fixed 
coordinate  system  [6],  by 


usHi  -  uSi  +  eijkQsjRsHk  +  dRsaJ  dt. 


(7) 


where  uSi  and  qSl  are,  respectively,  the  linear  velocities  along,  and  angular  velocities  about,  the 
Osxsixs2xs3  coordinate  system,  and  eiJk  is  the  permutation  symbol  defined  by 


^ ijk 


0,  if  any  two  of  i,  j,  k  are  the  same 
1,  if  ijk  is  an  even  permutation  of  1,  2,  3 
—  1,  if  ijk  is  an  odd  permutation  of  1,  2,  3. 


(8) 


In  standard  terminology,  <751  is  the  roll  rate,  qS2  is  the  pitch  rate,  and  qSi  is  the  yaw  rate.  The  relation¬ 
ship  between  the  angular  rates  and  rates  of  change  of  the  heading  angles  is  determined  [6]  from 

Qsi  =  Ru  dOsjldt,  (9) 

where  the  rotational  matrix  R  is 


R  = 


1 

0 

0 


0 

cos  05i 
-sin  05i 


—sin  052 
cos  052  sin  05i 

COS  052  COS  051 


(10) 
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The  transform  of  the  velocity  uSHi  from  the  ship  fixed  to  the  fluid-fixed  coordinate  system  gives 
the  velocity  Uom  in  the  latter  system  as 

u0 Hi~  TijUsHj.  (11) 

Applying  the  no-slip  condition  over  the  surface  of  the  ship  then  yields  the  boundary  condition  on  the 
fluid 

u0i~u0Hi  for  {*01,*02>*03)  =  {RoHl’RoH2’R()Hl}-  (12) 

In  addition,  assuming  infinitely  deep  and  unrestricted  waters,  we  have 

«o<  — '  «o~ /  as  V*o/*o/  —  oo,  (13) 

where  wo~/  *s  the  velocity  of  the  fluid  at  infinity  referenced  to  the  O0*oi*o2*03  coor<i*nate  system. 


Let  us  now  consider  conditions  at  the  free  surface  of  the  water.  We  designate  by  t)0(*oi,*o2T)  the 
elevation  of  this  surface  above  the  undisturbed  reference  plane  *03  =  0.  Two  conditions  prevail  at  the 
free  surface:  the  first  relating  to  continuity  of  stress  across  the  air-water  interface  and  the  second  relat¬ 
ing  to  continuity  of  the  interface  itself.  Neglecting  the  presence  of  any  ambient  surface  films,  the 
former  condition  can  be  written  [5]  as 


(n0  —  PA  —  2aHo)noj 


-d  -  rosV 


Bu0k 
dx0k  dx0i 


dUpy 

dx0j 


"ok  =  -vmnok 


at  *03  =  Vo- 


(14) 


Here,  PA  is  the  ambient  atmospheric  pressure  and  are  the  components  of  the  wind  stress  tensor  at 
the  surface.  The  coefficient  of  surface  tension  is  denoted  by  a.  Also,  «0/  and  H0  are,  respectively,  the 
outward  normal  to  the  free  surface  and  the  local  mean  curvature  of  the  free  surface.  The  latter  quan¬ 
tity  is  taken  as  positive  if  the  center  of  curvature  is  within  the  water.  From  Ref.  7,  we  have 


«o /  = 


9t)0 

9*oi 


8/1- 


9tjq 

9*02 


8/2  +  8,3 


1  + 

9t)0 

2 

j2 

9*oi 

+ 

9*02  J 

(15a) 


and 


2Hn  = 


92t?o 

9*oi 


1  + 


9t?o 


9*i 


02 


+  2 


9i7o  9t)0  d2vo 


9*oi  9*02  9*oi9*02 


92t?q 

9*o22 


1  + 


9t?0 


9*oi 


1  + 


9i7o 


Ox, 


01 


+ 


9tjo 


9*02 


3/2 


(15b) 


The  second  condition,  relating  to  continuity  of  the  air-water  interface,  is  obtained  as 
9t?0  ,  9tj0  ,97)o  A 

"9r  +  "M  9^T  +  “»!9^“"“_0  mxo>-V„-  (16) 

TRANSFORMATION  TO  MOVING  COORDINATES 

Suppose  that,  in  some  mean  sense,  the  origin  Os  of  the  ship  fixed  coordinate  system  is  moving  in 
the  negative  xoi  direction  with  velocity  U.  It  is  then  desirable  to  refer  the  governing  equations  of  fluid 
motion  to  a  coordinate  system  more  closely  aligned  with  the  position  of  the  ship  than  the  Oo*oi*02*03 
system.  We  define  this  system  as  the  Ox  1*2*3  system  with 

Xj  =  *0,  +  Ubat  (17) 
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and  we  have  taken  the  origins  of  the  two  systems  to  coincide  at  /  *  0.  The  relationships  between  the 
derivatives  in  the  fixed  and  moving  systems  are  found  as 


9  9 

9*0/  9*/ ' 


9 

9*i  * 


(18) 


We  also  have  that  any  quantity  Q0  in  the  fixed  system  is  given  in  the  moving  system  as 

Qo(x0ht)  =  0O(*,  -  U8nt,t)  =  Q(xjrt) 


(19a) 


and,  conversely, 


(?(*/,  0  =  Q(x  o/  +  U8nt,t)  =  Qo(x0i,t). 


(19b) 


If,  in  addition  to  the  above  transforms,  we  write  the  total  pressure  II  as  the  sum  of  its  hydrostatic 
and  dynamic  components  or 


n  =  pa  -  pgx 3  +  p, 

the  equations  of  fluid  motion  in  the  moving  coordinate  system  are  obtained  from  Eq.  (4)  as 

h (1 ~ +  u  a  ~ rs)  +  1(1 ~ ye)Uk]  =  0 

and 


fj  10  -  y «)«/)  +  U  7^-  {(1  -  yB)ui)  +  ((I  -  >«)«/»*} 


1  dP  ,  „  .  ,  1  9 

p  oXi  p  dxk 


0  -  rsV 


9m,  9m*. 


9x 


9x, 


2  9  m. 

T 8,;c  9T 


From  Eq.  (2),  the  void  fraction  yB  is 


ys~  J0  fvBVBdVB. 


(20) 

(21a) 


(21b) 

(22) 


The  gas  particle  velocity  uBi  and  location  RBi  =  R0Bi  +  U8nt  in  the  moving  coordinate  system 
are  determined  from  Eqs.  (3a)  and  (3b)  as 


and 


UBi  =  M,  +  URi  8,3 


(23a) 


—  RBi  —  U8n  +  uBi. 


(23b) 


To  resolve  the  boundary  condition  over  the  surface  of  the  ship,  it  is  convenient  to  first  define 

9Si  =  7tS,3  +  <t>Si  (24) 

so  that  the  ship-fixed  xS\  axis  is  basically  aligned  with  the  negative  xx  axis.  The  transformation  and 
rotational  matrices,  Eqs.  (6)  and  (10),  become 


—  cos  </>S3  cos  <t>s2  ~  cos  <f>S3  sin  <f>S2  sin  4>sl 

+  sin  <j>si  cos 


COS  05 3  Sin  052  COS  05i 
-  sin  053  sin  <f>S\ 


T  = 


-  sin  053  cos  052  -  sin  <f>S3  sin  052  sin  0Si 

—  COS  053  COS  051 


-  sin  05 3  sin  0S2  cos  05i 

+  COS  053  sin  051 


—  sin  052 


COS  052  sin  05! 


COS  052  COS  051 


(25a) 
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and 


Further,  we  write 


R  = 


0 

COS  051 

-sin  05i 


—sin  052 
cos  052  sin  051 
COS  052  COS  051 


uSi  =  -  UT,J  l6ji  +  \Si  =  -  UT,  1 1  +  vs/ 


(25b) 


(26) 


to  separate  the  basic  motion  of  the  ship  in  the  negative  xqi  direction  from  the  remaining  linear  velocity 
components  vSi  along  the  ship-fixed  axes.  The  velocity  of  a  point  on  the  surface  of  the  ship  then  is 
found  from  Eq.  (7)  as 

uSHi  =  ~UTaX  +  \Si  +  eiJkqSjRsHk  +  dR-sa/ dt-  (27) 

From  Eqs.  (11)  and  (12),  the  boundary  condition  over  the  surface  of  the  ship  with  respect  to  the 
Ox \x2x3  coordinate  system  is  obtained  as 


where,  from  Eq.  (5), 


ui  ~  TijUstjj  for  {x\,X2,xj}  —  {Rfi\,Rf[2,RHi} 


RHi  ~  RSi  +  TjjRsHj 


(28a) 

(28b) 


and  where  R$ ,  is  the  vector  from  O  to  0$.  In  addition,  we  have  from  Eqs.  (13)  and  (19a) 


^  — 1 -  uqoo  i  (x/  —  U8n(,t)  as  ^/X'X,  — 1 ►  00. 


(29) 


At  the  free  surface,  the  condition  on  continuity  of  stresses  is  found  from  Eqs.  (14)  and  (20)  as 


( P  ~  pgi)  -  2 aH)n,  -  (1  -  yB)/i 


Bu,  ,  Buk  2  »  Buj 

'  -  —"TO  a 


Bxk  Bx,  3  lk  Bxj 


(w) 

nk  =  “O' tk  nk  at  x3  =  rj 


(30a) 


where,  from  Eqs.  (15a)  and  (15b), 


8n  -  -r^-  8,2  +  8,i 


9xi 


Bx2 


1  + 


Or) 


9xi 


+ 


Or; 


u 


Bx2 


T7F 


and 


2H  = 


B2y 
Bx  12 

1  + 

Bx  2 

1  ^  dr)  drj  d2r)  d2r) 

Bx  1  Bx2  BxiBx2  Bx2 

1  + 

Brj  |2 
5xi 

1 

+ 

dr) 

2 

+ 

dr) 

213/2 

(30b) 


(30c) 
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The  condition  on  continuity  of  the  free  surface  is  obtained  from  Eqs.  (16)  and  (18)  as 

at  x3  =  rj . 


(31) 


Within  the  limits  of  the  inertialess  particle  hypothesis  and  the  neglect  of  ambient  stratification  and 
ambient  surface  films,  Eqs.  (21)  through  (31)  represent  an  exact  description  of  the  hydrodynamic  flow 
field  about  a  surface  ship.  The  set  of  equations  is  unsolvable. 

DEVELOPMENT  OF  THE  MODEL 


Part  1:  Decoupling  the  Fluid  and  Gas  Particle  Flow  Fields 


The  first  step  in  constructing  a  solvable  model  for  the  hydrodynamic  wake  of  a  surface  ship  is  to 
decouple  the  water  and  entrained  gas  particle  flow  fields.  This  entails  approximating  the  void  fraction 
by  yB  =  0  in  the  equations  of  motion  and  boundary  conditions.  Though  this  decoupling  step  is  not 
absolutely  essential,  it  does  greatly  simplify  the  numerical  work  involved  in  the  flow  field  calculations. 
If  higher  order  approximations  are  necessary,  an  iterative  scheme  such  as  outlined  in  Ref.  2  can  be 
used.  However,  since  the  equations  of  motion  in  the  inertialess  particle  limit  are  valid  only  for 
yB  «  the  error  introduced  by  setting  yB  =  0  is  generally  small.  This  is  especially  true  at  distances 
removed  from  strong  sources  of  gas  particles. 


With  yB  =  0,  Eqs.  (21a)  and  (21b)  reduce  to  the  standard  continuity  and  Navier-Stokes  equations 
for  an  incompressible  fluid,  or 


and 


(32a) 


|  jj  du,uk 

dt  9xi  dxk 


1  bJP  |  jit 

p  9x,  p  dxkBxk  ' 


(32b) 


Part  2:  Converting  the  Ship  and  Free  Surface  Boundary  Conditions  to  Known  Surfaces 

A  fundamental  difficulty  preventing  the  solution  of  Eqs.  (21)  through  (31)  is  that  the  ship  and 
free  surface  boundary  conditions  are  applied  at  surfaces  whose  locations  are  not  known  "a  priori." 
Hence,  the  second  step  in  constructing  a  solvable  model  for  the  hydrodynamic  wake  of  a  surface  ship  is 
to  convert  these  boundary  conditions  to  surfaces  whose  locations  are  known  "a  priori." 

a.  Ship  Boundary  Conditions 

Referring  to  Eqs.  (28b)  and  (25a),  we  see  that  both  RSh  the  vector  from  O  to  Os,  and  Tu,  the 
transformation  matrix,  contain  terms  dependent  on  the  instantaneous  location  of  the  ship  hull.  If  we 
stipulate  that  the  time-averaged  (mean)  motions  of  the  ship  about  its  equilibrium  motion  in  the  nega¬ 
tive  xqi  direction  are  zero  and  that  these  motions  are  themselves  suitably  small,  then  the  boundary  con¬ 
dition  on  the  ship  hull  can  be  applied  at  the  mean  location  of  the  hull.  For  those  points  on  the  ship 
hull  that  have  motions  in  addition  to  the  rigid  body  ship  motions,  the  boundary  condition  is  to  be 
applied  at  some  appropriate  mean  location  of  the  points  (for  example,  the  equivalent  actuator  disc  of  a 
propeller  [8]  or  the  zero  deflection  location  of  a  control  surface). 
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Taking  the  time  averages  of  Eqs.  (28b)  and  (25a)  and  imposing  the  above  stipulations,  we  find 
the  mean  location  of  the  ship  hull  as 

rhi  =  RSi  +  TjjRsHj,  (33a) 

where  overbars  denote  time  averages  (means),  and  where 

Tu  =  -8U  +  28/38y3.  (33b) 

To  the  same  order  of  approximation,  the  boundary  condition  at  the  mean  location  of  the  ship  hull  is, 
from  Eq.  (28a), 

Uj  =  Trusty  for  =  (34a) 

where,  from  Eq.  (27), 

uSHi  “  U&n  +  v5/  +  eiJkqSjRsHk  +  Vsrn-  (34b) 

Here,  dRSHi/ dt  =  v*SHi  indicates  that  effective,  as  opposed  to  actual,  velocities  may  be  required  at  those 
points  having  motions  in  addition  to  the  rigid  body  ship  motions.  (See  again  Ref.  8  for  a  discussion  of 
propeller  flow  problems.)  Combining  Eqs.  (34a)  and  (34b),  we  obtain 

ui  “  —Ubt\  +  Tjjivsj  +  ZjkiQskRsHi  +  vS//y)  f°r  {-^ i ^2^3}  =  (35) 


This  expression  is  equivalent  to  fhe  standard  linearized  ship  hull  boundary  condition.  Higher 
order  boundary  conditions,  based  on  Taylor  series  expansions  about  this  mean  condition,  have  been 
developed  [9].  However,  the  flow  fields  calculated  using  these  higher  order  boundary  conditions  do  not 
differ  significantly  from  those  calculated  using  the  mean  condition  while  the  numerical  work  involved 
does.  Hence,  we  take  Eq.  (35)  as  the  appropriate  ship  hull  boundary  condition  in  our  model  for  the 
hydrodynamic  wake  of  a  surface  ship. 

b.  Free  Surface  Boundary  Conditions 

The  exact  free  surface  boundary  conditions,  taken  at  the  unknown  surface  x3  —  77,  are  given  by 
Eqs.  (30)  and  (31).  If  we  stipulate  that  77  is  suitably  small,  these  boundary  condtions  can  be  approxi¬ 
mated  by  a  Taylor  series  expansion  about  *3  =  0.  To  model  the  hydrodynamic  wake  of  a  surface  ship, 
we  retain  only  the  lowest  order  terms  of  these  expansions— again  noting  that  higher  order  boundary 
conditions  lead  to  significant  increases  in  numerical  complexity  while  not  significantly  affecting  the  cal¬ 
culated  results. 


To  obtain  the  appropriate  boundary  conditions  to  be  applied  at  x3  =  0,  we  let  77  =  er 7'  where  e  is  a 
small  dimensionless  parameter.  Expanding  the  kinematic  free  surface  boundary  condition  given  by  Eq. 
(31)  about  x3  —  0  and  substituting  for  77,  wre  have 


e  +  U  +  u\  4-  677'  “t  “  +  0(e2) 
at  ox  3 

u3  +  *77'  —  —  +  0(e2) 

ox  3 


axj 


u2  T  677'  -7—^-  +  0(e2) 


dx-i 


dig' 

dx2 


=  0  at  X3  —  0. 


We  see  that  for  a  consistent  expansion  w3  =  ewj  and,  with  this  substitution,  find  the  kinematic  free  sur¬ 
face  boundary  condition  to  lowest  order  in  e  as 

+  (U  +  U\)  — 1-  u2  —  W3  =  0  at  x3  =  0.  (36) 

Ot  OXi  ox2 
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Consider  now  Eq.  (30a).  With  y#  =  0,  this  boundary  condition  becomes 
(P  ~  pgr)  —  2a//)  w,  —  <r iknk  =  —<rlknk  at  x3  =  17 

where  we  have  put 


v*  =  P 


Bui_  9»* 
dxk  9x, 


(37a) 


(37b) 


To  estimate  crik  near  the  free  surface,  we  note  that  from  continuity,  Eq.  (32a),  both  dui/dx\  and 
9w2/9x2  are  of  the  order  bu^/dx^  =  ebu^/bx^.  Hence,  cru  =  ccr/i,  <722  =  e(r22,  and  cr 33  =  ecr33.  Simi¬ 
larly,  since  there  is  no  reason  for  one  horizontal  direction  to  be  preferred  over  the  other,  cr12  =  ecr^. 
Also,  since  we  have  stipulated  small  surface  displacements,  both  bujbx3  and  bu2/bx3  must  be  small; 
or,  cr13  =  ecr{3  and  cr23  =  ecr23.  Thus,  we  find  cr,*  =  ecr'k.  Further,  from  Eq.  (30c),  we  have 


2H  =  —t 


9x,2 


9x,2 


+  0(63). 

Expanding  Eq.  (37a)  about  x3  =  0  and  substituting  the  above  results,  we  find 


:.P'  +  €2tj'  +  0(e3)  —  pg€7)'  4*  ae 

OX  3 


ay  !  ay 

+  0(€3) 

dx\  dx]; 

/  i  2 

« O’*  +  ez7) 


2  ,  9cr,* 


0x3 


+  0(c3) 


nk  =  -o-!kw)nk  at  x3  =  0 


where,  for  a  consistent  expansion,  we  have  set  P  =  eP'.  Taking  the  i  =  1,  2,  and  3  components  of  this 
expansion  and  observing  from  Eq.  (30b)  that  n\,  «2,  and  «3  are  of  orders  e,e,  and  1,  respectively,  we 
obtain,  to  lowest  order  in  e,  the  conditions  on  continuity  of  stress  across  the  free  surface  as 


P 


bui  9u3 


9x3  dxi 


P 


du2  9m3 


P  —  pgr)  +  a 


9x3  9x2 

2jl  a2 

dxi 


2  9x| 


=  o-!(3w)  at  x3  =  0, 

(38a) 

ate*  at  x3  =  0,  and 

(38b) 

„  a«3  ,  s 

2P  Bx  -  at  x3  —  0. 

(38c) 

Equations  (36)  and  (38)  are  the  desired  free  surface  boundary  conditions  for  our  model  of  the 
hydrodynamic  wake  of  a  surface  ship. 

Part  3:  Decomposing  the  Flow  Fields 

To  proceed  with  the  development  of  our  model,  we  next  wish  to  decompose  the  flow  field  into 
ship-induced  and  interaction-induced  components.  This  decomposition  allows  us  to  both  visualize  and 
calculate  the  overall  flow  field  as  two  basically  different,  but  interacting,  components. 

Let  us  write 


w,  =  v,-  +  v/  +  w,  (39) 

where  we  define  v,  and  v/,  respectively,  as  the  mean  and  randomly  fluctuating  flow  velocities  induced 
by  the  ship  motions  in  otherwise  still  water  and  w,  as  the  flow  velocities  induced  by  the  background  and 
background/ ship  flow  interactions.  Several  assumptions  are  implicit  in  the  decomposition  of  the  flow 
field  represented  by  Eq.  (39). 
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First,  we  have  assumed,  through  the  introduction  of  v/,  that  turbulence  modeling  will  be  incor¬ 
porated  into  our  description  of  the  hydrodynamic  wake  of  a  surface  ship. 

Second,  we  have  assumed  that  the  turbulent  velocity  fluctuations  v/  produced  by  the  ship  motions 
are  significantly  larger  in  magnitude  than  fluctuations  occurring  on  the  same  time  scale  in  the  back¬ 
ground.  Hence,  the  term  w\  has  been  neglected  in  the  velocity  decomposition  given  by  Eq.  (39).  This 
assumption  is  physically  reasonable  since  oceanic  sources  of  turbulence  are  normally  weak  compared  to 
the  sources  resulting  from  the  presence  of  a  ship. 

Third,  we  have  assumed  that  the  time  scale  of  the  turbulent  fluctuations  v/  is  small  compared  to 
the  time  scales  of  both  the  mean  ship-induced  flow  v,  and  the  mean  interaction  flow  wt.  This  assump¬ 
tion  permits  the  unambiguous  decomposition  of  the  velocity  field  given  by  Eq.  (39).  The  fundamental 
limit  on  this  assumption  is  the  time  scale  of  the  shortest  gravity-capillary  waves  present  in  the  back¬ 
ground.  For  example,  for  3.5-cm  waves,  we  find  that  the  turbulent  fluctuations  must  occur  with  a  fre¬ 
quency  >>20  Hz  which  is  near  the  limit  of  strict  validity  of  the  assumption.  For  much  shorter  waves, 
the  assumption  clearly  fails  though  certain  mitigating  factors  possibly  are  present  because  of  the  very 
small  depth  of  penetration  of  such  waves. 


In  a  manner  similar  to  the  decomposition  of  the  velocity  field,  we  write  the  dynamic  pressure  P 
and  surface  elevation  r\  as 

P  =  Pv  +  P'  +  Pw  (40a) 

and  .  f .  _  x 

7)  =  17v  +  TJv  +  Vwf  (40b) 

where  the  subscripts  v  and  w  refer,  respectively,  to  quantities  induced  by  the  ship  motions  in  otherwise 
still  water  and  to  quantities  induced  by  the  background  and  background/ship  flow  interactions.  The 
three  assumptions  applicable  to  the  decomposition  of  the  velocity  field  are  applicable  also  to  the  decom¬ 
positions  of  P  and  rp 

Proceeding  to  substitute  Eqs.  (39)  and  (40)  into  Eqs.  (32)  governing  the  fluid  motions  and  Eqs. 
(29),  (35),  (36),  and  (38)  governing  the  fluid  behavior  at  the  boundaries  of  the  flow,  and  time  averag¬ 
ing  the  resultant  expressions,  the  turbulent  equations  of  motion  and  the  turbulent  boundary  conditions 
are  obtained.  We  separate  these  equations  of  motion  and  boundary  conditions  into  sets  of  terms 
identified,  respectively,  with  the  flow  field  induced  by  ship  motions  in  otherwise  still  water  and  the  flow 
field  induced  by  the  background  and  background/ship  flow  interactions.  Setting  these  sets  of  terms 
individually  to  zero,  since  they  are  individually  separable  and  equal  to  zero  at  infinity,  produces  the 
governing  equations  of  motion  and  boundary  conditions  for  the  two  components  of  the  decomposed 
flow  field. 


a .  Ship-Induced  Flow  Field 


We  find  the  continuity  and  turbulent  Navier-Stokes  equations  for  the  ship-induced  flow  field  to  be 
given  by 


bxk 


=  0 


(41a) 


and 


9v,  9vf  9v,Vt 

~dr+  u  7^  + 

of  OX]  ox*. 


1  dPy 


+ 


p  dxj  dxk 

The  boundary  conditions  applicable  to  this  flow  field  are  obtained  as 


3v,  — 

- — - V.Vju 

P  dxk 


Ml. 


V| 


0 


as  v  x,x,- 


(41b) 

(41c) 
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and 


P,  ~  PSV\  +  « 
drj 


9vj  9v3 
9x3  9x! 

9v2  9v3 

— -  H - - 

9x3  9x2 


92i?v  +  92T) 


9x,2 


9x2 


+  vW 

for 

{xi,X2,X3}  — 

(4  Id) 

+  Or'3v  =  0 

at  x3  —  0, 

(41e) 

+  o-23v  =  0 

at  x3  =  0, 

(4  If) 

_ 

dv3 

< 

+  crT'v  —  2  jji 

dx3 

—  o-33v  =  0  at  x3  =  0, 

(41g) 

+  U  ~z~~  +  (v,Tjv  +  v/i^)  —  v3  =  0  at  x3  =  0. 


(41  h) 


dt  '  dx\  '  9x, 

In  developing  Eqs.  (41e)  through  (41g),  we  have  allowed  for  the  possibility  that  the  average  values  of 

both  the  fluctuating  stresses  a'ijv  and  fluctuating  surface  tension  <rT'v  are  not  zero.  Also,  Eq.  (41h)  has 
been  developed  by  writing  Eq.  (36)  in  the  equivalent  form 

*L+UgL  +  !ZL-Ul.0 

ot  dx\  9x, 

which  follows  from  Eq.  (32a)  and  the  fact  that  17  is  independent  of  x3. 
b.  Interaction-Induced  Flow  Field 


For  the  interaction-induced  flow  field,  the  governing  equations  of  motion  and  boundary  conditions 
are  found  as 


and 


9w,  dw: 

T7+u^~L  + 
ot  9x] 


Pw  -  PSVw  +  a 


dVw  +  u  dy* 


dt 


9x] 


dwk 

T±  =  0’ 
dxk 

(42a) 

d  ,  .  1  dPw  n  92w>( 

(w,wk  +  WjVk  +  v.Wk)  _  +  g-  r  ,  , 

QXk  p  ox,  p  dxkdxk 

(42b) 

■*  /(•*/  —  Uhnt,t)  as  V x,x,  — » 

(42c) 

0  for  {xi,x2,x3}  =  {PhI’Ph2>Ph^ > 

(42d) 

dwi  dw3 

dx3  dxi 

=  o-1(3k,)  at  x3  =  0, 

(42e) 

8  w2  9  w3 
9x3  9x2 

n 

q 

09 

5 

n 

0 

(42f) 

92tjh 
9xj2  9x2 

“j  2F  “  -o'3(3m’)  at  x3  =  0, 

(42g) 

wiVv  +  v,T -  w3  =  0  at  x3  =  0. 

(42  h) 

Part  4:  A  Matched  Asymptotic  Expansion  for  the  Ship-Induced  Flow  Field 

Equations  (41)  describing  the  ship  induced  flow  field  in  otherwise  still  water  are,  with  currently 
available  techniques,  unsolvable  on  large  domains  regardless  of  the  turbulence  model  used.  The  princi¬ 
pal  difficulties  preventing  a  solution  are  the  nonlinear  free  surface  boundary  condition  given  by  Eq. 
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(41  h)  and  the  second  order  viscous  derivative  terms  in  Eq.  (41b).  These  derivative  terms,  which  are 
influential  only  in  the  vicinity  of  the  ship  hull  and  the  ship  turbulent  wake,  give  Eqs.  (41)  the  charac¬ 
teristics  of  what  Van  Dyke  [10]  terms  a  singular  perturbation  problem.  For  such  problems,  further 
progress  can  be  made  by  applying  the  method  of  matched  asymptotic  expansions  [10]  to  the  solution 
procedure.  Before  applying  this  method,  we  must  nondimensionalize  Eqs.  (41)  to  identify  the  relevant 
governing  parameters,  and  before  this,  introduce  the  rudiments  of  a  turbulence  model. 


For  our  present  purposes,  it  is  sufficient  to  use  isotropic  eddy  viscosity  and  diffusivity  models  for 
the  turbulent  stresses  and  diffusivities  appearing  in  Eqs.  (41).  More  complex  models  can  be  readily 
incorporated  without  changing  the  basic  results  of  this  section.  Following  Rodi  [11],  we  write 


~T~i  Q'ikv 

-V/V*  =  — 


9v,  dvk 
'  +  * 


9x, 


dx, 


—  y  ks ik> 


(43a) 


where  v,  is  the  turbulent  viscosity,  and  k  is  the  turbulent  kinetic  energy  per  unit  mass, 


k  =  y  v/v/. 


By  analogy  with  Eq.  (43a),  we  also  take 


<rTV  ®  / 


92t?v  +  9lv 


9*1  Bx{ 


(43b) 


(44a) 


where  a,  denotes  the  coefficient  of  turbulent  surface  tension.  Further,  in  analogy  with  eddy  diffusivity 
models  for  turbulent  heat  or  mass  transport,  we  put 

-v/tk  =  o-vv,  (44b) 


where  (tv  is  a  dimensionless  number.  We  also  set 


and 


T"  "  v  ~ 

P 

a  =  (Taat. 


(45a) 
(45  b) 


Here,  v  is  the  kinematic  viscosity  of  the  water;  crv  is  the  ratio  between  this  and  the  turbulent  viscosity; 
and  cra  is  the  ratio  between  the  coefficient  of  surface  tension  and  the  coefficient  of  turbulent  surface 
tension. 


Since  our  primary  interest  is  in  the  hydrodynamic  wake  produced  by  the  ship,  we  choose  nondi- 
mensionalizing  parameters  associated  with  this  wake.  As  a  length  scale,  we  use  the  dominant  Kelvin 
wavelength  [12]  27rt/2/gand,  as  the  time  scale,  the  dominant  Kelvin  frequency  lirU/g.  This  choice  of 
scaling  parameters  renders  the  wake  Froude  number  equal  to  unity.  The  dimensional  quantities  appear¬ 
ing  in  the  ship  wake  problem  then  are  given,  in  terms  of  their  dimensionless  counterparts,  as 


2nU 

t  —  -  T, 


2  it  U2  , 

*/  =  — —  £i> 


V/  =  Uvit 


pv  =  p  U2PV,  vs,  =  UvSh  v*sw  =  Uv 


SHi’ 


Qsi  — 


2t tU 


5  R  Mil  5 

QSi>  KSHi  ~  KSHi> 

o 


n  -  27 tJj1  n 

KHi  —  K 


Hi> 


p  _  2t tUI  - 
Rsj  R  Si>  and  rjv 

6 


2t t UI 

g 


(46a) 
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For  the  turbulent  kinetic  energy,  we  take 


k  =  ^k 
2t tU 


(46b) 


which  insures  that  both  dimensionless  components  of  Eq.  (43a)  are  of  the  same  order  of  magnitude,  as 
required  for  dimensionless  analysis. 

Substituting  Eqs.  (43)  through  (46)  into  Eqs.  (41),  we  obtain  the  dimensionless  equations  govern¬ 
ing  the  ship-induced  flow  field  problem  as 


*-0. 

96 


(47a) 


9  6  96_  dv,vk 
9r  96  96 


BPV 


+ 


96  9£ 


R 


Kt 


av/  3v,  Bvk  2  - 

<T>,  TZ - r  — - r  ” - T  kbik 


v  96 


96 


96 


(47b) 


v,  —  0  as  V66 


(47c) 


v,  =  —8,i  +  Tjj(vSj  +  eJkiqSkRsHi  +  Vsty)  f°r  {6>6>6l  =  {R  h\>R  hi>R  Hit  <  (47d) 


9vj  9v3 


96  96 


=  0  at  6  =  0. 


(47e) 


9v2  9v3 


96  96 


=  0  at  6  =  0. 


(47f) 


Pv  —  2nf)v  +  ( aa  +  1) 


aKt 


R 


Kt 


and 


92rjv  +  92tjv 


962 


96 


2te,  +  1)  85,  +  2  J_  _  „  at{j_0  (4?g) 


R 


Kt 


96 


3  R 


Kt 


9rjv  9^  9 


dr 


96  96 


V/TJv 


CTt,  9rjv 


96 


V3  =  o  at  6  =  0. 


(47h) 


In  these  expressions,  Rja  and  aKl  are,  respectively,  the  turbulent-flow  Reynolds  number  and 
turbulent-flow  inverse  capillary  number,  both  based  on  the  Kelvin  wake.  We  have 

»  _  27TC/3 

KK!  ~ 


gVt 

OLt 


277  t/3  j. 

<rv  - =  o’ VR^ 


aKt  ~ 


gv 

a 


O’  V 


ock 


(48a) 

„  n  (48b) 

pUvt  c ra  pup  aa 

where  RK  and  aK  are,  respectively,  the  laminar-flow  Reynolds  number  and  laminar-flow  inverse 
capillary  number.  Typical  values  of  RK  and  aK  for  water  are 

Rk  =  6.4  x  10 5U\  aK  =  73/  U 

with  U  in  m/s  (1  m/s  —  2  knots). 
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As  explained  in  Ref.  10,  the  method  of  matched  asymptotic  expansions  is  a  solution  technique  for 
problems  such  as  given  by  Eqs.  (47)  where  the  higher  order  derivatives  have  only  limited  influence  on 
the  flow  field.  The  technique  consists  of  developing  two  series  solutions,  one  governed  by  an  outer 
flow  parameter  that  takes  into  account  the  unimportance  of  the  higher  derivatives  for  this  flow;  the 
other  governed  by  an  inner  parameter  that  recognizes  the  limited  influence  of  these  higher  derivatives. 
The  final  or  composite  solution  is  obtained  by  an  appropriate  combination  of  the  two  series  that  is  valid 
over  the  entire  flow  field. 

The  terms  outer  and  inner  have  arisen  historically  because  of  the  application  of  the  method  to 
boundary  layer  and  wake  flows  in  infinite  domains.  For  these  problems,  it  is  easy  to  visualize  two  dis¬ 
tinct  flow  regions.  For  the  ship  wake  problem,  as  we  shall  see,  this  analogy  falters,  because  both  outer 
and  inner  solutions  occupy  the  flow  region  aft  of  the  ship.  However,  the  formalism  of  the  method  does 
not  depend  on  defining  specific  flow  regions,  and  we  retain  the  terms  outer  and  inner  only  for  historical 
consistency. 


a.  Ship-Induced  Outer  Flow  Field 


Following  the  formalism  laid  out  in  Ref.  10,  we  seek  an  outer  asymptotic  solution  to  Eqs.  (47)  for 
large  RK(=RKl/o-v)  since  the  inverse  of  this  parameter  multiplies  the  highest  order  derivatives  in  the 
equations.  We  write 


v,(r, £/',/?£ )  ~  [  lim  8CR*)]  v/(r 

J 

(49a) 

Py( T.frM  ~  1  lim  8(^)1  PHt.O. 

(49b) 

ic(r, £,;/?*)  —  f  lim  8 (/?*)]  F(t,£,), 

(49c) 

K—  J 

>£\>£i\Rk)  ~~  [  Am  SCR*)] 

(49d) 

where  a  superscript  t  denotes  the  flow  variables  associated  with  the  outer  flow  field,  and  where  8  (Rx) 
is  an  expansion  parameter  to  be  determined.  Substituting  these  expansions  into  Eqs.  (47)  and  taking 
the  limits  as  Rx  — *  we  find 


dv/ 

a?7 


-o, 


dv/  dv/ 

dr  e*. 


f  lim  S(RX) 


aw 


&P! 


(50a) 

(50b) 


v 


t 

i 


lim 


1 

8(RX) 


dv/  dv/ 

e?7+  877 

dv/  dv/ 
df 3  d^2 

Py  —  2irvjy 


0 

0 

0 


v/—  0  as  —  oo, 

f  5/1  +  Tfj(vsj  +  ejkldSk^SHI  +  vj#/)] 

for  **  [R  H\>R  Hl’R  Hi)' 

at  f  3  -  0, 

at  ?3  -  0, 


(50c) 


(50d) 

(50e) 


at  £3  -  0, 
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and 


dvt  dj?v 

9t 


lim  8(Rk) 
rk~°° 


—  v/  “  0 


at  £3  =  0. 


(50h) 


We  see  from  Eq.  (50b)  that  the  highest  order  derivatives  have  been  lost  in  the  outer  asymptotic 
expansion  as  has  kf.  The  outer  problem  is  thus  reduced  to  an  inviscid  flow  problem  (we  set  kf  =  0 
since  it  is  indeterminate).  Having  lost  the  highest  order  derivatives,  the  boundary  conditions  on  the 
ship  hull  given  by  Eq.  (50d)  cannot  be  identically  satisfied  nor  can  the  entirety  of  free  surface  stress 
conditions  given  by  Eqs.  (50e)  through  (50g).  We  instead  specify  zero  fluid  penetration  through  the 
hull  or,  equivalently,  zero  normal  velocity  and  continuity  of  the  normal  stress  component  across  the 
free  surface.  Hence,  Eqs.  (50e)  and  (50f)  are  ignored  in  the  outer  flow  approximation.  Additionally, 
we  drop  Vsw  from  the  outer  problem  since  these  effective  velocities  are  associated  with  vorticity  pro¬ 
duction,  which  is  inconsistent  with  the  inviscid  nature  of  the  outer  problem. 


Defining  the  normal  to  the  ship  hull  with  respect  to  the  ship  fixed  coordinate  system  by  nSHi ,  the 
normal  nt  in  the  Ox  1X2X3  coordinate  system  becomes,  consistent  with  the  suitably  small  restrictions 
already  imposed  on  the  ship  motions, 

ni  =  TijflsHj •  (51) 

The  appropriate  outer  flow  boundary  condition  on  the  ship  hull  then  is  obtained  as 


= 


lim 


8C Rk) 


[  n\  +  n,  Tjj  {vSj  +  eJklqSkRsHi )] 


for  {fi,£2>£3}  —  {Rm>R Hi)- 


(52) 


From  this  equation,  we  note  that  the  only  nondegenerate  outer  flow  problem  results  from  assigning  a 
finite,  nonzero  value  to  lim  8(RK).  Without  loss  of  generality,  we  can  set  this  limit  to  unity. 
rk~°° 


The  outer  flow  problem,  since  it  is  inviscid  and  irrotational,  can  be  simplified  considerably  by 
introducing  a  dimensionless  outer  flow  velocity  potential  f  related  to  its  dimensional  counterpart  (f>t 
through 


Then 


<f>f  = 


2tt  £/3 

g 


4>f. 


(53a) 


MI 

'  ’ 

and  Eq.  (50a)  yields  Laplace’s  equation  for  the  potential  as 


KMt 


=  o, 


(53b) 


(54a) 


while  Eq.  (50b)  can  be  integrated  to  give  Bernoulli’s  equation 

Py  +  +  \{  +  \  v/v/  =  0.  (54b) 

or  z 

The  appropriate  boundary  conditions,  which  all  can  be  cast  in  terms  of  <j>f,  are  given  by  Eqs.  (50c), 
(50g),  (50h),  and  (52). 


The  above  problem  is,  essentially,  still  unsolvable  because  of  the  nonlinearities  appearing  in  Eqs. 
(50h)  and  (54b).  However,  as  pointed  out  by  Newman  [12],  neglecting  these  nonlinearities  is  necessary 
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for  consistency  with  the  approximations  already  made  in  reducing  the  free  surface  boundary  condition 
to  the  plane  f3  =  0.  This  is  easy  to  verify.  Recall  that  for  this  reduction,  we  took  r}/  ~  ei }/'  where  e 
was  a  small  dimensionless  parameter.  This  gave  v/  ~  ev/'  and  PJ  ~  eP$'  in  the  vicinity  of  the  free 
surface.  The  free  surface  boundary  condition  was  then  expanded  about  |3  =  0,  and  all  terms  of  higher 
order  than  e  were  neglected.  Referring  to  Eq.  (54b),  we  see  also  that,  for  the  outer  flow  problem, 
—  ev/'  in  the  vicinity  of  the  free  surface.  Hence,  4>t  and,  consequently,  v/  are,  respectively,  of 
order  and  ev/'  in  the  vicinity  of  the  free  surface.  Thus,  the  nonlinear  terms  in  Eqs.  (50h)  and 
(54b)  are  of  order  e2  for  the  outer  flow  problem  and  must  be  dropped  to  ensure  consistency  with  our 
previous  approximations. 


Returning  to  dimensional  quantities  via  Eqs.  (46a)  and  (53a),  we  can  summarize  the  ship-induced 
outer  flow  field  problem  as 


eV 

dxfixj 


=  0, 


v/ 


_ 

as  -JXiXi 


(55a) 

(55b) 

(55c) 


and 


v/n,-  =  -Unx  +  «/ Tjj (ySj  +  ejklqSkRSH,)  for  {xi,x2,x3}  =  {RH{,RH1}RH3), 

+  2U -f}— +  U2  +  g  T^-  =  °  at  x3  =  0, 

Bt 2  BtBxx  dxf  Bx3  3 


together  with  the  prognostic  relations 


Rj  = 


Bt1 


Bt 


+  U 


B<j/ 


Bx  i 


(55d) 

(55e) 


(55f) 


and 


,T=- 


8 


&_  +  u&_ 

Bt  9x! 


at  x3  =  0. 


(55g) 


Equations  (55)  are  recognized  as  being  the  traditional  Kelvin-Neumann  problem  [12]  for  determining  a 
ship’s  Kelvin  and  radiated  wave  systems.  This  problem,  while  nontrivial,  is  solvable. 


One  cautionary  note  is  in  order.  For  the  outer  problem  to  be  fully  consistent  over  the  entire  flow 
domain,  it  is  necessary  that  the  n\  component  of  the  ship  normal  be  everywhere  small  (of  order  e). 
For,  if  at  some  points  on  the  hull  it  is  of  order  unity,  we  find  from  Eq.  (52)  that,  at  those  points, 
v/  ~  0(1)  as  opposed  to  order  e.  Hence,  the  consistency  of  the  outer  flow  approximation  breaks  down 
at  such  singular  points.  Physically,  such  points  give  rise  to  the  nonlinear  bow  wave  problem  [12]  which 
must,  in  itself,  be  treated  by  other  expansions  and  then  incorporated  within  the  traditional  Kelvin- 
Neumann  problem. 


b.  Ship-Induced  Inner  Flow  Field 

Having  lost  the  viscous  derivatives  in  the  outer  asymptotic  expansion,  we  have  been  forced  to 
neglect  ship  boundary  layer  and  propeller  effects  and  the  ship  turbulent  wake  generated  by  these  effects. 
These  viscous  phenomena,  which  are  confined  to  thin  regions  adjacent  to  the  ship  hull  and  circum¬ 
ferential  to  a  line  (the  resultant  thrust  axis)  running  aft  of  the  ship,  must  be  recouped  in  the  inner 
asymptotic  solution.  The  technique  leading  to  this  solution  uses  the  thin  nature  of  the  viscous  regions, 
and  the  concomitant  knowledge  that  flow  quantities  vary  much  more  rapidly  across  these  regions  than 
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along  other  regions,  to  rescale  and  simplify  the  equations  and  boundary  conditions  governing  the  gen¬ 
eral  ship-induced  flow  field. 

Since  our  interest  is  in  the  hydrodynamic  wake  of  the  ship,  we  begin  by  considering  the  inner 
asymptotic  solution  for  the  region  aft  of  the  ship.  Before  starting,  however,  one  more  cautionary  note 
is  in  order.  The  subsequent  developments  restrict  the  resultant  thrust  axis  to  be  parallel  to  the  direction  of 
motion  of  the  ship,  that  is,  canted  neither  significantly  towards  nor  away  from  the  free  surface. 

Following  again  the  formalism  laid  out  in  Ref.  10,  we  introduce  an  inner  expansion  parameter 
A  (Rk)  that  represents,  in  dimensionless  terms,  the  circumferential  extent  of  the  viscous  wake  region 
about  the  resultant  thrust  axis.  This  parameter  requires  the  property 

A  (Rk)  *  0  as  R %  *  °° .  (56) 

Since  the  region  is  thin,  a  nonuniform  scaling  of  coordinates  is  necessary  to  ensure  that  derivatives 
across  and  along  the  region  are  a  uniform  length  scale.  We  accomplish  this  by  modifying  the  dimen¬ 
sionless  global  coordinates  f  /  as 

=  s\,  £2  =  &(RK)s2,  £3  —  AC Rk^s3  (57) 

where  the  sf  give  the  required  uniform  inner  coordinates.  We  have,  here,  tacitly  assumed  that  the  wake 
axis  parallels  £1,  from  which  assumption  the  above-noted  restriction  on  the  direction  of  the  resultant 
thrust  axis  has  been  obtained. 


Based  on  the  scaled  coordinates,  we  seek  an  inner  asymptotic  solution  to  Eqs.  (47).  By  definition, 
the  dimensionless  velocity  along  the  wake  axis  must  exist  if  a  long,  thin,  viscous  region  running  aft 
of  the  ship  is  to  exist.  Hence,  must  be  relatively  independent  of  the  circumferential  extent  of  the 
wake,  and  we  write 

~  vf(r,s,)  (58a) 

where  a  superscript  *  denotes  the  flow  variables  associated  with  the  inner  flow  field.  The  velocities  v2 
and  v3  transverse  to  the  wake  axis  are  taken,  more  generally,  in  terms  of  the  inner  expansion  parameter 
as 


and 


v2(t, $;/?*) 


lim  A x(Rk) 

rk~‘°° 


v3(t,Si',Rk) 


lim  A k(RK) 

rk-°° 


where  the  nonnegative  exponent  \  is  to  be  determined. 


(58b) 

(58c) 


Substituting  Eqs.  (57)  and  (58)  into  Eq.  (47a),  we  obtain  the  continuity  equation  for  the  inner 
flow  field  as 


3v  * 
dsi 


+ 


lim  Ax~HRk) 


9v*  9vf 

+ 


ds2  9s3 


=  0. 


(59) 


Thus,  for  the  fundamental  property  of  conservation  of  mass  to  be  preserved  in  the  inner  flow  field,  we 
must  take  X  =  1.  Employing  this  result  and  Eqs.  (57)  and  (58)  in  the  dimensionless  version  of  Eq. 
(43a),  we  then  find 

k(r, $,•;/?*)  ~  k*(T,s,)  (60) 

if  both  components  of  the  inner  form  of  Eq.  (43a)  are  to  be  of  the  same  order  of  magnitude,  as 
required  for  dimensionless  analysis. 
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Using  Eq.  (48a),  Eqs.  (47b)  become 


9v  * 

9r  9$!  9s, 


9vf  9vfv, 

+  ^  .  =  lim 


a^v  ,  9 

ds\  ds\ 


cruR 


(c Tv  +  2) 


9vf 
9  Si 


1 
3  K 


9S2 


1 


<r„/?*A20?*) 


,  9v*  ,  9v$ 

(o-„  +  1)  1 -  +  A20?*) 


9s2 


9si 


+ 


9s3 


1 


*prkaHrk) 


9vf  9v? 


(61a) 


9vf  9v*  9v*v* 

-r- 1  +  T^  +  “^=  lim 
9r  9  Si  os,  rk~°° 


a^v  9 


A20?*)  9s2  9si 


1 


V  vRk^2(Rk) 


4=(^)(„  +  !,«£!  +  9£i 


9si  9s2 


+ 


9s2 


1 


i>Rk&2(Rk) 


(  avi  2  -» 

+  2)  a77  “  T k 


+ 


9  s3 


9v|  9v3  9v3v? 

+ -r2-1  =  lim 
9r  OS]  9  s, 


<tvrkaHrk) 


1 


+,)«!+  a*5 


9s3  9s2 


_ a^v  _j_  9 

A2(RK)  9s3  9s! 


ctvRkA2(Rk ) 


A2(R/')(fTv  +  1)  .  + 


(61b) 


9v* 


9si  9s3 


+ 


9s2 


o'  irRjsr^O?*) 


\c, 


9v3  9v* 


+  1)  ^-z-  + 


9  s2  9  s3 


+ 


9s3 


1 


<tvrkaHrk) 


9v?  9  - 

^  +  2)-^f-fk' 


(61c) 


Examining  these  expressions  in  the  limit  /?*  — 1 ►  °° ,  we  see  that,  for  the  inner  asymptotic  expansion  to 

be  both  well  behaved  and  to  retain  viscous  flow  effects,  lim  RKA2(RK)  must  be  a  finite,  nonzero 

rk^°° 

number.  Hence,  without  loss  of  generality,  we  set 


which  yields 


lim  RkA2(Rk)  =  1 
rk~°° 


AO?*)  =  \/Rp. 


(62a) 

(62b) 
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With  this  limiting  value  for  RKh2(RK),  we  see,  in  addition,  that,  to  eliminate  singular  terms  from  the 
inner  asymptotic  equations, 


lim  A 2(Rk) 


P*(T,Sj). 


(63) 


Using  the  identities 


and 


9 

1  9v? 

_  Bv*  q 

i 

9  s, 

o-„  9s2 

9  s2  9  s, 

<r„ 

1  9v? 

Bv*  a 

1 

— 

9  s, 

crv  9s3 

9$3  dSj 

<*V 

(64a) 


(64b) 


which  follow  from  Eq.  (59)  and  assuming  (soon  to  be  demonstrated)  that  the  wake  region  is  steady  in 
the  frame  of  reference  moving  with  the  ship,  the  inner  flow  equations  for  describing  this  region  reduce 
to 


9v* 
9  s,- 


=  0, 


(65a) 


9v*  9vf v* 

+ 


9sj 


9  s, 


B 

Bs2 


9v*  9v2v* 


9si 


+ 


9  s, 


9s2 


+ 


cr„  +  1  9v* 

+  j_ 

<r„  +  1  9vf 

cr„  9s2 

9s3 

<rv  9s3 

;  +  Qvr  9 

>  0s2  dSi 

1 

vv 

9 

cr„  +  1  9v2 

1  il 

<tv  +  1  9v2 

9s2 

cr„  9s2 

3  <jv 

+  — 
9s3 

O',  9^3 

(65b) 


(65c) 


9si 


3  +  9v|v* 


9  s, 


9P;  9v*  Q 

9s3  9s3  9s, 


O’,, 


+ 


9 

crv  +  1  9v* 

o-v  +  1  9v| 

* 

<N 

9s2 

crv  Bs2 

9s3 

C T„  9s3 

3  crv 

(65d) 


We  note  that  all  second  order  derivatives  involving  sb  along  with  axial  gradients  of  the  pressure  and 
turbulent  kinetic  energy,  have  disappeared  from  the  inner  flow  problem.  Hence,  Eqs.  (65)  are  hyper¬ 
bolic  in  the  Sj  direction  and  elliptic  in  the  (s2,s3)  cross  plane  (i.e.,  parabolic  overall).  Mathematically, 
this  means  that  the  values  of  the  inner  flow  variables  on  some  plane  s\  +  ds\  are  determined  entirely 
from  their  values  on  the  plane  s\  together  with  the  boundary  conditions  on  the  plane  s\  +  ds\.  Physi¬ 
cally,  this  implies  that  if  the  values  of  the  inner  flow  variables  are  steady  on  some  plane  Si  aft  of  the 
ship,  they  remain  steady  as  we  progress  further  downstream  from  the  ship.  Since  we  have  already  re¬ 
stricted  the  unsteady  motions  of  the  ship  about  its  mean  motion  to  be  small,  the  primary  sources  of  the 
wake  are  the  boundary  layer  of  the  ship  in  uniform  motion  and  the  commensurate  propeller  effects. 
Both  of  these  sources  are  nominally  steady;  thus,  our  assumption  that  the  wake  region  is  steady  in  the 
frame  of  reference  moving  with  the  ship  has  been  validated. 
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Let  us  turn  our  attention  now  to  the  boundary  conditions  for  the  inner  flow  problem.  We  assume 
that  somehow  (computationally,  experimentally,  or  parametrically)  the  values  of  the  inner  flow  vari¬ 
ables  on  a  plane  aft  of  the  ship  have  been  specified.  Hence,  as  follows  from  the  parabolic  construction 
of  the  inner  flow  problem,  the  boundary  conditions  on  the  ship,  given  by  Eqs.  (47d)  are  irrelevant 
regarding  the  subsequent  downstream  values  of  these  variables.  Substituting  Eqs.  (57),  (58),  and  (63) 
into  the  remaining  boundary  conditions  [Eqs.  (47c)  and  (47e)  through  (47h)],  using  Eqs.  (48)  and 
(62),  and  invoking  the  steady  nature  of  the  wake  region,  we  find 
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In  arriving  at  the  latter  two  boundary  conditions,  we  have  taken 
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where  the  nonnegative  exponent  k  is  to  be  determined. 
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Consider  Eq.  (66d) .  For  k  less  than,  equal  to,  and  greater  than  unity,  respectively,  we  have  in 
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(68c) 


The  condition  obtained  for  k  <  1  is  independent  of  the  inner  flow  field  and,  except  for  the  solution 
tjv  =  0,  gives  a  continuously  increasing  elevation  for  s2  >  0  since  [(cra  +  lVcrJa*  is  positive. 
Hence,  we  determine  either  k  >1  one  <  1  with  y  J  =  0.  With  either  constraint,  Eq.  (66e)  yields,  in 
the  limit  RK  —*  °°, 


v$  =  0  at  S3  =  0.  (69) 

Together  with  this  result,  Eq.  (68c)  overspecifies  the  boundary  conditions  on  the  inner  flow  problem. 
Thus,  we  require  for  a  well-stated  problem  either  k  =  1  or  k  <  1  with  y*  =  0.  Since  the  former  con¬ 
straint  is  less  restrictive,  we  conclude  k  —  1. 
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Summarizing,  the  boundary  conditions  for  the  inner  flow  problem  are  ascertained  to  be 

v*— 0  as  V s22  +  S32  -*°°  (70a) 

and 

3v*  3v* 

v*  sTT “  "aTT “ 0  ats’-°  <70b) 

along  with  the  prognostic  relation  for  ij*  given  by  Eq.  (68b).  We  note  from  Eq.  (70b)  that  the  free 
surface  acts  as  a  plane  of  symmetry  for  the  inner  flow  field.  Hence,  we  arrive  at  the  important  result 
that  the  solution  to  the  ship-induced  inner  flow  problem  is  identical  to  the  lower  half  space  (s3  <  0)  solution 
that  would  be  obtained  by  solving  the  inner  flow  problem  for  the  ship  and  its  image  (about  Sj  =  0)  in  an 
infinite  fluid. 


Returning  to  uniform  dimensionless  quantities  via  Eqs.  (57),  (58),  (60),  (63),  and  (67)  and, 
thence,  to  dimensional  quantities  via  Eqs.  (46),  we  obtain,  with  the  aid  of  Eqs.  (45)  and  (48),  the 
dimensional  form  of  the  ship-induced  inner  flow  field  problem  as 
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c.  Ship-Induced  Composite  Flow  Field 

The  ship-induced  composite  flow  field  is  derived  most  straightforwardly  by  the  method  of  additive 
composition  detailed  in  Ref.  10.  Basically,  to  obtain  the  uniform,  first  order,  composite  solution,  we 
sum  the  first  order  outer  and  first  order  inner  solutions  and  subtract  the  part  they  have  in  common  so 
that  it  is  not  counted  twice.  Formally,  this  procedure  translates  into  summing  the  first  order  outer  and 
inner  solutions  and  subtracting  either  the  outer  expansion  of  the  inner  solution  or  the  inner  expansion 
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of  the  outer  solution.  Because  we  have  developed  a  matched  asymptotic  solution,  these  latter  two 
expressions  are,  by  construction,  equal  and,  in  our  case,  identically  zero. 

To  see  this,  let  us  rewrite  the  inner  problem  in  terms  of  the  independent  outer  variables  Sub¬ 
stituting  Eqs.  (57)  into  Eqs.  (65)  and  (70)  and  taking  the  limit  as  RK  —*  °°,  we  find  the  outer  represen¬ 
tation  of  the  inner  flow  problem  as 
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where  the  superscript  *  t  denotes  the  outer  value  of  the  inner  variable.  Hence,  the  outer  expansion  of 
the  inner  solution— that  is,  the  solution  of  Eqs.  (72)— is 


v?(g,)  =  0. 


(73) 


For  completeness,  let  us  also  rewrite  the  outer  problem  in  terms  of  the  independent  inner  variables  st. 
Substituting  Eq.  (57)  into  the  nondimensional  versions  of  Eqs.  (55)  and  taking  the  limit  as  RK  — 1 ►  °°, 
we  find  the  inner  representation  of  the  outer  flow  problem  as 
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where  the  superscript  t  *  denotes  the  inner  value  of  the  outer  variable.  Here,  the  boundary  condition 
on  the  ship  hull  has  been  neglected  since  the  matched  asymptotic  solution  encompasses  fully  only  the 
region  aft  of  the  ship.  Thus,  the  inner  expansion  of  the  outer  solution  satisfies  a  two-dimensional 
Laplace’s  equation  with  zero  gradient  boundary  conditions.  From  Ref.  13  and  Eq.  (74b),  we  obtain 
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Consequently,  the  method  of  additive  composition,  together  with  Eqs.  (58),  yield  the  ship- 
induced  composite  flow  velocities  aft  of  the  ship  as 
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or,  in  dimensional  terms,  as 
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Let  us  now  turn  our  attention  to  the  ship-induced  composite  pressure  and  surface  elevation  fields. 
From  Eqs.  (63)  and  (67),  both  the  inner  pressure  and  elevation  are  of  second  order  in  the  inner 
expansion  parameter  h(RK).  Hence,  they  do  not  contribute  to  the  uniform,  first  order,  composite 
pressure  and  elevation  fields.  However,  the  inner  flow  velocities  do  have  a  first  order  corrective  effect 
on  the  outer  elevation  field  and,  as  a  result,  on  the  outer  pressure  field. 


To  demonstrate  this,  we  return  to  the  dimensionless  statement  of  the  kinematic,  free  surface 
boundary  condition  for  the  ship-induced  flow  field  given  by  Eq.  (47h).  Substituting  Eqs.  (76)  into  the 
kinematic  condition,  dropping  the  products  of  v/  and  t)v  since  they  are  of  order  e2  in  the  surface  eleva¬ 
tion,  and  using  Eq.  (70b),  we  find,  in  the  limit  RK  —  °o, 
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This  expression  represents  the  uniform,  first  order,  kinematic  free  surface  boundary  condition. 
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where  £v  gives  the  correction  to  the  outer  elevation  field,  we  obtain,  with  the  aid  of  Eq.  (50h),  the 
governing  relation  for  £v  as 
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This  correction  to  the  outer  dynamic  head  leads  to  a  proportionate  correction  in  the  outer  pressure 
field.  We  have 
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which,  since  PJ/vJ  is  well  behaved,  is  itself  well  behaved. 
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Returning  to  dimensional  variables  via  Eqs.  (46),  we  find 

TJv  =  1?/  +  C v 


and 
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SUMMARY 

A  model  for  the  hydrodynamic  wake  of  a  surface  ship  has  been  developed.  Neglecting  only 
ambient  stratification  and  ambient  surface  films,  we  have  begun  with  the  general  equations  and  bound¬ 
ary  conditions  for  the  flow  field  and,  through  a  series  of  rational  approximations,  reduced  the  problem 
to  a  solvable  one. 


23 


R.  A.  SKOP 


The  approximations  are: 

•  The  fluid  dynamics  are  independent  of  the  gas  particle  (bubble)  dynamics. 

•  The  ship  undergoes  only  small  motions  about  a  uniform  mean  motion. 

•  The  free  surface  elevations  are  small. 

•  The  ship-induced  turbulence  is  significantly  stronger  than  the  ambient  turbulence. 

With  these  approximations,  the  overall  problem  becomes  separable  into  one  governing  the  ship-induced 
flow  in  otherwise  still  water  and  one  governing  the  interaction  of  the  background  with  the  ship-induced 
flow. 


The  former  problem  was  examined  further  by  applying  the  method  of  matched  asymptotic  expan¬ 
sions.  We  found  that  the  ship  induces  an  outer  flow  field  that  satisfies  the  traditional  Kelvin-Neumann 
problem  for  calculating  a  ship’s  Kelvin  and  radiated  wave  systems.  This  outer  flow  field  is  determined 
by  Eqs.  (55).  We  found  also  that  the  ship  induces  an  inner  flow  field— its  turbulent  wake.  This  flow 
field  is  governed  by  the  steady,  three-dimensional,  parabolic  Navier-Stokes  equations  with  plane  of 
symmetry  boundary  conditions  at  the  free  surface.  The  inner  problem  is  given  by  Eqs.  (71).  The 
overall,  or  composite,  ship-induced  flow  field  is  derived  from  these  two  asymptotic  flow  fields  via  Eqs. 
(77)  and  (82). 

Besides  the  plane  of  symmetry  boundary  conditions  at  the  free  surface  for  the  turbulent  wake 
problem,  the  more  significant  results  of  this  paper  include: 

•  Surface  tension  effects  are  unimportant  in  determining  the  uniform,  first  order,  ship- 
induced  flow  field  in  otherwise  still  water;  and 

•  The  inner  flow  field  produces  a  first  order  modification  to  the  traditionally  calculated  ship 
wave  elevations.  This  result  could  explain  the  experimentally  observed  differences 
between  the  transverse  wave  systems  of  model  and  full-scale  ships. 
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